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Thiemobwectiememortethis Study 1s to develop classification 
Sugeest onsc tO £LoOnecast the daily probability of occurrence of 
marine fog at selected locations on the west coast of the 
United States, using parameters easily obtainable from upper- 
ae SCcunagings and Ssuriace observations. In order to achieve 
this objective a computerized stepwise linear discriminant 
analysis program is extensively employed. Data input con- 
Sists of surface and radiosonde observations for the five- 
vedr pen toqmmmumn~ lt 765 LO SO June 1973 at three U. S. west 
coast stations, namely San Diego and Oakland, California and 
Surilayute, Washington. 

Tables showing the number of fog and no-fog cases, the 
elasS@i@earmron Lunctions, and the percentages of correct fog 
and no-fog discrimination are presented for each station. 
The most capable fog/no-fog discrimination parameters are 
discussed for each set of classification equations. Test 
results for the San Diego equations using a three-year inde- 


pendent data set are also shown. 
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ie INE RODUCTION 


Marine fog is an impediment to both commercial and mili- 
Pate comm Oo pind wcnGdealt re Operations. Such fog has often caused 
costly delays to commercial ship owners and domestic air 
carriers while endangering lives and property as a result of 
owe Sibi litvtes. Naval @perations such aS aircraft transfer 
and landing, antisubmarine warfare and amphibious maneuvers, 
navigational reconnaissance and search and rescue missions 
may be hampered by marine fog. The importance and impact of 
marine fog on naval operations during World War II and the 
recent costs (1969-1974) in lives and property damage of 
ship and aircraft accidents suffered by the United States 
Navy during fog situations has been described by Wheeler 
(O74)... 

Recognized aS an important meteorological problem, marine 
fog prediction and dissipation have been the subjects of on- 
going research by several Department of Defense (DOD) acti- 
vities. In particular the Naval Air Systems Command (NAVAIR) 
and the Office of Naval Research (ONR) have jointly sponsored 
marine fog research in an attempt to model and simulate the 
POC@muCGrapUrDOSsSes Of Drediction. A Naval Postgradute School 
Ce comoeamp under the direction of Dr. Leipper, Chairman “of 
Prem ilewaeement Ol Oceanography, and Dr. Renard of the Depart= 
Memon eMeteOorealogy, 1S participating in this coordinated 


Girone with the following specific objectives: 
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Wie LOmobserve and describe the formation and dissipa- 
tion TOL fog at sea in cooperation with Cornell 
Aeronautical Laboratory, Inc. 

2. to relate the synoptic variables in meteorology 
and oceanography which are important in marine fog 
formation and dissipation to microphysical and 
mMicrometeorological characteristics. 

So.) Tesimprove marine fog forecasting and to develop 
marine fog prediction aetinedis. "= 

Several recent NPS studies have addressed the problem of 

MoIwines tog Climatology Of the North Pacific Ocean and incor- 
porated a method of synthesizing ships' surface synoptic re- 
ports for the purpose of deriving frequencies of marine fog 
Stwanenecm amg tebretson, 1974; Daughenbaugh, 19/75; Renara,; 
pigeememetson and Daughenbaugh, 1975; Wilims, 1975). Wallace 
(1975) and Hale (1975) have explored the feasibility of uSing 
meteorological satellite data as a means of discerning the 
presence of marine fog. Pod stuaywoy Mee lure (1974) aoa spo; 
tential link between temperature inversions and the occur- 
rence of marine advection fog was investigated, in order that 
temperature inversion parameters could possibly be used as 
forecasting aids at coastal stations. 

Petnougie thie Environmental Prediction Research Facility 


(EPRF), Monterey, California, is engaged in numerical 


Presta arey beet wane others, 1973: Observation, Analysis 
and Prediction of Marine Fog. Naval Postgraduate School, 
Department of Oceanography Annual Report for period 1 July 
ieaneo 50 syune 1973),, p. 10. 
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mode beng Studies on the physics of fog formation, mainte=- 
manee "ana Gissipation (Feit, 1972; Barker, 1973), and 
researchers elsewhere, such as Mack et al (1973), have 
aeveloped numerical models of advection fog, few opera- 
tional forecasting schemes for marine fog on the west 
coast of the United States have been developed to date. 
Fleet Numerical Weather Central (FNWC), Monterey, Cali- 
frmumba, Nas an Operational advection fog program ("“FTER") 
which forecasts grid point fog probabilities to 48 hours 
based on a system of multi-parameter tests (U. S. Naval 
Weather Service Command, 1975). Given the correct humid- 
ity values, this program gives acceptable results in cases 
Ob Dbroad-scale fog. However, itS primary weakness stems 
from the large 63x63 grid -- local effects just cannot be 
reconciled. 

"Forecasting fog is one of meteorology's hardest prob- 
lems. This is because there are so many factors involved-- 
moisture, temperature, wind, geography, stability, cloud 
cover, and several mieten OH Le yrs hoped thatethnis study 
may provide insight and methodology which will aid fore- 
Casters tn predicting the occurrence Of marine fog fEor 


Seuteec rea locations Om the west coast of the United States: 


25. S. Naval Weather Service Command, 1975: Wien S . 


Naval Weather Service Numerical Environmental Products 
Manual, NAVAIR , Department of the Navy, 
Wastindeon, §D. C., p. 4.26=2. 
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18 Oi BACKGROUND 


A. eR yee WwW OFMSERECTED PREVIOUS FOG RESEARCH 

RmcOnscmactanle amoumt of research into the problem of 
marine fog on the North American west coast has been con- 
@uetred in the southern California area. In 1948 Leipper 
(1948) devised a model to illustrate the manner in which 
Hog Stiuations GQevelop ian the San Diego area. This model 
was described in four stages and is still used to some de- 
Gecemovyemilitary £orecasters in the southern California 
area. 

jmhe first Stage or initial conditions in Leipper's model 
requires the presence of air over the sea which has a tem- 
perature higher than the sea-surface temperature and which 
is relatively dry aloft. This condition guarantees the for- 
MattoneoLean inversion which will restrict the vertical 
movement of moisture and thus cauSe a thin lower layer of 
Tearoom approach Saturation. Such initial conditions are 
typically brought about when a lobe of the North Pacific sub- 
Pmeptical anticyclone pushes inland over northern California 
causing a general easterly flow over southern California so 
Mismneive walt arriving at San Diego is warm and dry, Maving 
deccenqea Ehe 4000-foot slope of the coastal mountains. 

In the second stage of fog development, the easterly 


flow weakens, the flow of warm subsiding air toward the sea 
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Peemeases, anc the air offshore from San Diego remains in a 
Mematively sStagnant condition. Due to evaporation from the 
Bea the moisture content of the lower layers increases con- 
Siderably, and the heat content decreases due to the conduc- 
Peomeor heat downward. During this process a surface 
inversion is formed and the lowest air layers become nearly 
Poimaarcead at a temperature close to that of the underlying 
Sea Surface. 

Fog forms in the third stage as the normal northwesterly 
airflow and sea breeze regime returns. The fog is created 
in the thin surface layer at sea as the warm, nearly satu- 
rated air moves over the colder coastal waters. LMS ye fos 
Pomorocucea by cooling from below (Leipper, 1968). Once the 
Weiee Cloud is formed, incoming radiation reflected from the 
Popmeana Meat radiated from the cloud itself results in cool- 
ing of the thin fog layer to a temperature lower than the 
GSeeam Stisrtace in the vicinity. Then as the sea breeze 
strengthens during the daytime, the fog is advected shore- 
ward. 

Pres ftourth and final stage 1s characterized by the deepen 
ing of the nearly adiabatic lower layer to approximately 
400 feet or more. Since evaporation from the sea continues 
and the presence of the strong inversion still restricts 
EmMcmevereErcal Movement Of moisture, the dew point remains high 
in the marine layer. After the daytime sea breeze carries 
pHemucGgmuover the shore, radiation and conduction of heat £rom 


the land mass dissipate the fog somewhat, but the nearly 
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meager c layer, being deeper than the previous stage, is 
Hememaestroyed as rapidly. Rather, it retains its high 
moisture content in vapor form and moves inland. As the 
moemne layer beneath the inversion continues to deepen, a 
point is reached where the mixing and cooling processes do 
Heme result im condensation throughout the layers Then the 
Mpper portion Fills with water particles first and fog 
gees not form on the ground. In this way the fog sequence 
is ended and the stratus regime begins. The complete fog 
sequence Lelipper described usually extends over a period of 
apeut five days. 

In a later work Leipper (1968) compared the meteorologi- 
cal conditions associated with the observance of a sharp 
smog bank near Riverside, California with those previously 
shown to be related to winter fogs in southern California. 
He believed that the characteristics of certain stages of 
these fog situations in winter were quite similar to those 
Secervedupy es saringer (1963) in July 1961 and to those re- 
Poea@eedmapyeotephens (1965). The distinguishing characterzs-— 
tic was the presence of unusually warm, dry air over the 
inversion and avery cold marine layer beneath it. In such 
Situations the air above the inversion is often nearly 10 C 
warmer than the underlying sea. The combination of very 
Wem air aloft and an unusually cold marine layer results 
Mmeehcestrongly Stable situation described by both Edinger 
and Stephens as having high smog potential when it moves in- 
land. Such situations in the San Diego area were associated 


Wait ogs havang sharp boundaries. 
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Three nondiurnal indices presented by Leipper in 1948, 
aie cOounad useful in the prediction of west coast fog, also 
Hayepe Important in predicting situations favorable to the 
Shallow sharp-banked smogs which are often observed. The 
indices are composed from data obtained in the morning radio- 
sondes used together with the coastal water temperature and 
the surface dew point in the marine layer. Leipper (1948) 
described the indices as follows: 

le Height of the Inversion BaSe: The height above 
which the air temperature increases with height at 
me MOsSt rapia rate on the morning raob (radivosonde 
eoSservatilon) . 

2. Temperature Index: The quantity (Ta-Tw); if an in- 
version exists with base below 3,000 feet, Ta is the 
highest air temperature above the inversion base on 
ee Morning raob, otherwise Ta 1s the suriace air 
temperature, and Tw is the coastal sea~surface tem- 
PebaQeuke. 

BS. e tiemeture Index: The difference between the afternoon 
dew-point and the coastal water temperatures. 

Himcmeconarttons £Ound most favorable for Leg occurrence, 
Tf applicable to smog, require the base of the inversion to 
bemobelkew sl SOO Eeet;, the temperature index to be positive, 
Cnemenemltorstiure index to be greater than —-5 C. 

Gehtocder et al (1967) describe three significant aspects 
Sueieaine alr invasion along the west coast of the United 


Srotes, two Of which are related to fog development. The 
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most noticeable invasion is that of the diurnal sea breeze. 
This well known phenomenon is produced by differential heat- 
mice or the land and water masses and results in the sea 
beeeze circulation cell which brings cool, moist marine air 
Momawara. One of the most distinctive characteristics of 
miomoca oreeze 1S the £front-like appearance at its leading 
edge. 

The inland movement of the sea breeze and its front is 
Maomoly controlled by terrain features. In addition, the 
large-scale weather patterns affect the sea breeze. When the 
Evm@ptiLe-scale gradient flow is directed offshore, the sea— 
breeze front is intense but does not penetrate far inland, 
and the strength of the flow is weak. However, when the 
gradient flow is onshore, the front is weak and the sea- 
breeze flow is strong, resulting in considerable inland 
penetration. 

Perhaps the most significant marine air invasion is that 
of the United States west coast monsoon. This phenomenon be- 
Gens e2neethne late spring and continues until fall, bringing a 
slow, steady transport of marine air inland from the North 
Pete een Ocean high pressure cell. This subtropical high 1s 
composed of a shallow marine layer capped by a subsidence 
inversion. The marine inversion is the most pronounced and 
shallow during the monsoon season. trregulacs £inectuatione 
in the marine inversion are generally related to the synoptic 
weather patterns, while the diurnal variations of the inver- 
wHoOmeaepth are partly related to the interaction of the sea 


breeze and the monsoon. The monsoon is generally confined to 
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pmeeeet teCSS than 2 km deep and its inland path is primarily 
through low gaps in the coastal mountains. 

In summary, the monsoon is a feature of the general 
Manos mmherle: Clrembation which undergoes modification at the 
coastline, interacts with both the sea breeze and the moving 
synoptic-scale systems, and is confined to a very shallow 
layer. 

One of the most comprehensive studies of marine-fog fore- 
Gasting on the west coast was done at Scripps Institution 
of Oceanography during a contract with ONR, and was super- 
vised by D. F. Leipper (Leipper and others, 1948). The 
purpose of this "Fog Project" was to develop and to test 
principles which might serve as a basis for fog forecasting 
in coastal areas and possibly elsewhere. 

The "Final Report of the Fog Project" outlines some of 
he theory of condensation and evaporation which have parti-= 
cular importance in the development and application of the 
fog forecasting technique. It is generally believed that 
Enrom locals tog Lorecasting problem is primarily one of pre- 
Sectimgmthe moisture content of the air and the amount of 
Pcmeootleng which will occur. The sea-surface temperature 
is an important variable in determining the amount of modi- 
fication (measured by the change in content of heat and 
moisture) which will be brought about in a given overlying 
air mass. This temperature is also a convenient quantity 
with which other more erratic gquantities such as air tempera- 
ture and dew point may be compared, as was done in Leipper's 


Gloete ye nen-diurnal fog indices. In low coastal regions the 
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dew-point temperature is a reliable index of moisture con- 
tent, and is the most convenient quantity for use in fore- 
fasting changes in moisture which will affect fog formation, 

The forecasting procedure devised by the researchers of 
the "Fog Project" consists of three parts. First the fore- 
Casterpe muStuprepare pregnostic “charts im order to determine 
from them certain fundamental weather features (local fea- 
tures such as wind speed, cloud cover, range of ceiling and 
surface pressure; also, synoptic-scale features such as 
changes in surface air flow and subsidence). The second step 
Mmoives applying statistical aids (e.g., climatological data, 
temperature and time forecasts, and visibility data) to deter- 
Mine the local changes in heat and moisture content of the 
air associated with the fundamental features which are fore- 
cast. Finally, the forecaster must determine the visibility 
from the moisture content, changes in heat content and other 
available data by applying empirical rules. 

Imebormation obtained from radiosondes such as the height 
Of the anversion base, the strength and thickness of the in- 
Wemclom metic. , May be significant in fog forecasting, in the 
Jinaiweregect” §. Bs Wickham discusses why the inversion in= 
dexmeneitance Of the Inversion base) iS important. Generally 
speaking, when a strong temperature inversion exists over a 
nearly adiabatic layer, the moisture flux through the inver- 
Broneic somailnas Compared with the flux through the surface 
layer. Consequently, moisture is trapped below the inversion 
amamover a period of time, near saturation of the entire 


Baudbpatic bayer Can occur. Not only does near saturation 
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Seemt most rapidly with low inversions, but the amount of 
cooling necessary to produce condensation throughout the 
Weaewet 1S less for thin layers than for thick ones. Thus, 
ene time required to approach saturation and the amount of 
Boolmng needed to form fog are both functions of the height 
Seeehe Inversion base. The inversion index is of further 
importance Since it effects the diurnal temperature range 
and indicates the location of the dry marine~-air boundary. 

A more recent investigation into the relationships be- 
tween temperature inversions and the occurrence of marine 
advection fog was undertaken by McClure (1974). He studied 
Meee primary Inversion parameters: the height of the in- 
version base, thickness of the inversion layer, and the tem- 
perature gradient within the inversion layer, using ten 
months of surface and upper-air data (July 1973 to April 
1974). The upper-air sounding stations used were San Diego 
Hiemeqometyeriteld), California; Oakland (International Air- 
poate), sCautrOornia; and OQuillayute, Washington. The surface 
Observation stations were San Diego (Lindberg Field), Cal- 
iowa sOoakhand (International Airport), California; and 
Seattle-Tacoma (International Airport), Washington. 

Because Of the considerable geographic separation 
(nearly 100 miles) between the upper-air soundings taken at 
Oumibavyiitece and the surface fog observed at the Seattle-Tacoma 
Report, no définitive relationships were achieved. However, 
at San Diego a relationship between the change in inversion 


Parameters and the occurrence of fog was found to exist in 
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Been che summer and winter seasons. At Oakland, consistency 


in the change of One parameter was shown for most of the 


Gata period. 


A brief summary of some of McClure's relationships between 


temperature inversions and marine advection fog follows. 


ile. 


Summer fog cases are generally preceded by an inver- 
Sion in the atmosphere below 850 mb. For soundings 
taken duptnguerog. the lack of an inversion below 

850 mb is the exception rather than the rule, 
especially during the summer and early fall. 

Most fall and winter fog cases (80%-90%) are pre- 
ceded by an inversion. 

During the summer season in San Diego the thickness 
of the 0000 Greenwich Mean Time (GMT) inversion 
layer usually decreases within 24 hours prior to the 
occurrence of fog and the inversion base lowers. 
Cooling of the surface layer precedes the fog. Where 
fog persists there is a general warming trend in the 
lower atmosphere associated with a subsidence inver- 
Sronwewin the OOOO GMT sounding. 

During the winter season in San Diego the inversion 
base rises in the last 24-hour change of the 0000 
GMT sounding before a fog Sequence. in the -L2Z003Gi 
sounding the inversion is based at the surface for 
Celeste memour soundings prior to the fog. Within 
PoenourseOoL: the commencement of fog, the strength of 
thewinverston layer bessens, the depth increases, 


while the temperature gradient decreases. 


on 





Ss In Oakland's summer Season there is a decrease in 
the inversion thickness and an increase in the tenm- 
perature gradient in the 0000 GMT sounding within 
24 hours prior to the fog occurrence. 

6. In Oakland's fall season there is almost always a 
surface inversion in the last two 1200 GMT soundings 
Pweloumec Ene Onset Of fog. This 1s also true for 
the winter months of January and February. 

1 During the winter season in Oakland inversions do 
not alwayS accompany or precede fog. If an inversion 
Sersoes, it iseusually weakening prior to fog occur= 
rence, while the inversion thickness and temperature 
gradient are decreasing. 

McClure's investigations and recommendations Served as 
the platform from which this study was launched. A much 
broader data base and computer analysis was believed to be 
essential if more concrete relationships were to be made. 
Additional parameters were needed, especially within the 
surface observations. Finally, it was considered to be 
highly desirable to have the upper-air sounding station and 
the surface observation station at the same location. These 
changes and others, plus a different analysis approach, were 


MoecOorpOonated in this research. 
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tie OBbVECTIVE, DATA AND DEFINITIONS 


A. OB@GECTIVE AND APPROACH 

The objective of this study was to develop classifica- 
Mmromeedidatltons to forecast the daily probability of occeur— 
rence of marine fog at selected locations on the west coast 
of the United States, using parameters easily obtainable 
from upper-air soundings and surface observations. In order 
to achieve this objective, a stepwise linear discriminant 
analysis program, BMDO7M (Dixon, 1973), and the IBM 360 com- 
puter were extensively employed. Application of these power- 
ful computational resources to the available data would 
hopefully statistically reveal what parameters are most 


Srandticant in marine-fog prediction. 


B. Date DESCRIPTION 

Upper-air soundings and surface 3-hourly airways reports 
for the five-year period of 1 July 1968 to 30 June 1973 were 
used for three U. S. west coast stations: San Diego, Cali- 
Ponta Oaktand, California; and Quillayute, Washington, 
Phieselnatial data base included almost 11,000 soundings (twe 
Sa monaeaver 43,700 ‘suriace reports (eight daily}. Sea- 
Setace temperature information for the period of interest 
Wascmextracted from the Fishing Information Bulletin Supple-— 


Memes published by the U. S. Department of Commerce. The 


data provided by the Naval Weather Service Detachment, 
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Asheville, North Carolina were taken from the National Cli- 
ftoeerc Center*s historical data files, known as Tape Data 
Family-14 and -56 (TDF-14 and TDF-56), Airways Surface Obser- 


vations, and Radiosonde Data. 


C. DEE INETION OF FOG AND S2@G DAY 

For purposes of this study "fog" (as in marine advection 
fog) is defined as a visible aggregate of minute particles 
of water (droplets) based at the earth's surface, which re- 
feecemior at Zontalevyicibl VityecOo less titan seven miles and te 
Smemstruction to vision at the reporting station. A "fog 


day" as used in this study was defined as the period from 
fooler st on one Calendar day through 1600 PST on the follow= 
Imemecalendar day, so that a day would usually include both 


Micmtormatilon and QGissipation time of fog (after Leipper, 


19). 


24 





JU = EXPERIMENTAL PROCEDURES 


A. DATA PREPARATION AND ORGANIZATION 

Pn anyestudy anveolving large amounts of information per— 
maps Che most difficult and crucial first step is the proper 
Peemeratlon and Organization of the data. Initially the 
data tapes had to be screened for missing reports and modi- 
Filed for the use of FORTRAN programming on the NPS IBM 360 
Gomputer, All surface and upper-air reports were identified 
Byeeyeas, Month, day and hour numeric time groups. The upper-— 
alr sounding time group identifiers were converted from 
Greenwich Mean Time (GMT) to Pacific Standard Time (PST) in 
order to conform with the airways surface observation time 
Groups. The sea-surface temperature (SST) data were extrac- 
EeCmnvEeema month (lst to the 15th and 16th to the 30th/2ist) 
PeoOnMmimeomosctechakrLts in’ the Fishing Information Bulletin Sup= 
Prememecm Uo. Department of Commerce, 1968-1973), and were 
gagdeq tontne surface report data. The temperature informa- 
PEOUNOMetmhe SUbface reports was then converted from the 
PaiGbennelt scale to the Celsius scale, thereby utilizing the 
MeoLimewsystem fOr all data fields. 

Bacmwes-hnourly surface observation contains 35 data Eields 
(not including the tape and station number and the date), 
Viechwmdesctribe the following meteorological parameters: 


CoMbbtniggeyescibijilty, wind direction and speed, dry-bulb 
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temperature, wet-bulb temperature, dew-point temperature, 
melatiuve humidity, sea-level pressure, station pressure, sky 
=omartlon, Cloud information (amount, type, and height of as 
Many as four cloud layers plus the total amount and total 
opaque amount, as well as summation amounts at the second 
and third layers), and other atmospheric phenomena, includ- 
Pin iiid Storms, liquid and frozen precipitation, and “6b- 
Semuletrons tO vision. This last data field is especially 
—migoimtamt: Since the definition of a fog occurrence is based 
upon the presence of fog at the station as an obstruction 
Peer orou we rurthermore, the code lists three types of Log: 
ogee sce £09 wueand “qround fog". Although marine advec- 
MiOnieoge is not explicitly indicated, the great majority of 
mogmocecurrences were described as "fog" and are probably 
marine fogs considering the three stations selected for this 
Steud Vy « 

The upper-air soundings were reported twice daily (0400 
and 1600 PST) and contained the following data for both man- 
datory and significant levelS: pressure in millibars (mb), 
height in meters (m), temperature in degrees Celsius, rela- 
eivemmumidilty, and wind direction and speed. Since fog is 
a surface phenomenon, the entire sounding was not needed. 
Momeeporecdupny Mececlure (1974), the important inversion in= 
homnmaeionmrelating to the occurrence of fog existed at 850 mb 
and below. Consequently, each upper-air report was arbi- 
pecietivEruncatled at 700 mb. 

igemoutiace Observations and the upper-air reports were 


integrated on one tape by numeric time groups and organized 
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Peyeebocataon. The five years of data were subdivided into two 
Seasons, designated "dry" and "wet". This classification was 
prompted by the fact that west coast weather, especially 
California, has two dominant seasons characterized by a 
generally dry summer period and a wet winter period. The 
so-called “dry" season consisted of the six months from April 
through September while the remaining six months, October 
through March, comprised the "wet" season. Initial attempts 
to study the data by months instead of by seasons was unsuc- 
cessful because the monthly sample size of fog cases was 
Sften too small. (This was especially true for the summer 
months at Oakland, California.) 

Nineteen meteorological parameters, derived from both 
the surface and the upper-air data, were selected for analy- 
sis. Since upper-air observations were made only twice a day, 
the nineteen parameters were calculated for the times of the 
upper~air soundings, namely 0400 PST and 1600 PST. Never 
Perec ee Sour rye suridce report was checked for Ehe oceur- 
rence of fog and such information was retained in the 
poeateter ew tOGgC strength’, FS. (A discussion of the parameters 
mor OunGde im esSsection B.) 

As the parameters were calculated, the reports were fur- 
Pieimemlassiited Into Lour categories, based’ on the occurrence 
Of Mom—oceuLrence Of fog over the previous three days, and 
whether fog was reported within the current 24-hour period. 
Pmm—eacgqumurreo fog) occurred between 1601 PST on one calendar 


fveancde GOO PST on the following calendar day, that period 
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Toeeelassiuredmas aw fog” day ("no-fog" day). If fog was 
Mopermecanon any Surface Observation in the 72 hours (3 days) 
peamgretO the Gurprent ELog/no-fog day, that period was termed 
MeOgeniStOruywe. rE, On the other hand, no £fo0g was reported 
cmereng the bast 72 hours, a "“no-fog history" period had oc= 
Suerecde se USiNGsehese definitions, the four categories shown 
in Table I were derived. Generally speaking, categories 1 
and 2 depict the possible onset of fog since they have no- 
fog histories, while categories 3 and 4 depict the possible 
Geeopersion Of £LOg Since they have fog histories. 

The surface observations and upper-air data were now in 
a form which could be easily analyzed. Application of the 
discriminant analysis technique required the data to appear 
in two or more groups. The categories with common histories 
and opposite current fog designations were concatonated and 
PeocoenneoebMeoyM im order to develop discriminant functions. 
For example, categories 1 and 2 (3 and 4) were used together 
Pere Groum 1(3) consists of “foggers” and group 2(4) con- 
SrSts One no-Loggers" in the current 24-hour period. 

ie baca Glmitatitons and Special Considerations 

In any set of data there are always certain problems 

Oce limitations which must be reconciled or at least recognized. 
Mivemmest ObVEOuS Kestriction in this study was that the data 
covered only a five-year period (July 1968 to June 1973) and 
MmeomGieniurer just three west coast stations. These stations 
eremsecparated by hundreds of miles and extend over 15 degrees 


Otearatrtude. Also, recall that the surface observations were 
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mace at three-hour intervals (0100, 0400, 0700, 1000, 1300, 
feo, 1900 and 2200 PST), while the upper-air soundings were 
feumetyred Only twice a day (0400 and 1600 PST). It is impor- 
Pameteto realize that although fog may have occurred during 
mee hOUr Period; Say at 1000 PST, a full list of the para= 
fase Sea fOr EOgG discrimination was available only at the 
peumeaing times, O400 and 1600 PST. Considerable amounts of 
fog were reported at times intermediate to the upper-air 
Tarr nis implies that significant changes in the existing 
sounding, especially in the marine layer, May Neve occured 
immer cwr interim, resulting in the formation of fog. These 
subtle changes would not be reflected necesSarily in either 
Mime mp per—alr Or surface observations for that fog day. 

Several other limitations of the data set should be 
iemenowmed.  htrst Of all, if any information field in either 
the surface observations or the radiosonde soundings was 
missing in a reporting period (0400 or 1600 PST), that re- 
Pomseeiadecowbe discarded. Unfortunately, missing data or a 
Missing report caused a gap in the 72 hours of "fog history” 
and resulted ina loss of three days of information. Never- 
theless, only about five percent of all the reports were 
rejected as a result of missing data. Another restriction 
which the data set imposed upon this research was that the 
BLOdgadeSYNOpElc picture was not considered. That 1s, sounding 
Stations were not used to derive advective quantities or 
define air-mass boundaries, nor were synoptic charts 


employed. Furthermore, the sea-surface temperature (SST) 
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Gataetaken from the Fishing Information Bulletin Supplements 
may not always be the most representative SST nearest the 
pwEeace Station. Although the SST changes rather slowly, a 
Dbi-weekly average SST can be inaccurate, especially in the 
Summer months when solar radiation can result in significant 
temperature changes. The Oakland station SST was a special 
problem because the station is located on the east side of 
San Francisco Bay, Several miles from the ocean regime. The 
sea-surface temperatures nearest the bay outlet were used, 
Since the marine fog which forms offshore and affects Oak- 
land generally enters through this opening. 

One final special consideration was recognized in 
San Diego where the surface observation station (Lindberg 
Field) and the upper-air station (Montgomery Field) are not 
colocated. The radiosonde launch site is about 5.4 miles 
inland from Lindberg Field and is elevated 124 meters above 
Sea Mevel (or about 115 meters above the surface station). 
Under these circumstances a low-level inversion within the 
124-meter layer could exist and yet not be shown on the 
BouncdumcgemeMeclure felt that in fog situations where inver=— 
Sions did not exist, but seemingly should have, there may 
have been an inversion near the sea surface. Consequently, 
in this study a special procedure was implemented to “bogus 
Miroir ace Tnversion under the proper circumstances. Le 
the surface temperature at Montgomery Field (upper-air 
peat wonvmewasmwequal to OF greater than the surface tempera— 


mime wateLandberg Freld (surface observation station) an 1so- 


thermal condition or low-level inversion was said to have 
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exrsted. fi V~thne sounding had no inversion shown, the thick- 
mess Of the low-level inversion was taken to be the differ- 
ence in elevation between Montgomery and Lindberg Fields 
erz4a4m-9m=ll5m). It is recognized that this is an artificial- 
fey and that Other natural influences may account for these 
differences, such as the nature of the underlying surface 

and the difference between the marine air and urban environ- 


ments. 


Ee DISCUSSION OF THE FOG/NO-FOG DISCRIMINATION PARAMETERS 

The selection of the meteorological parameters to be 
used as fog discrimination variables was facilitated by pre- 
vious researchers in this subject. The three non-diurnal 
indices developed by Leipper (1948) as fog forecasting aids 
Piesoanm Diego were adopted with only minor modifications for 
mse an this study. McClure’s thesis (1974) provided informa- 
tion about three important inversion parameters: height of 
the inversion base (also one of Leipper's indices), thickness 
of the inversion layer, and the temperature gradient within 
the inversion layer. 

UWsineGmene wOrks of Leipper and McClure as guidance, nine=- 
teen initially selected parameters were chosen or computed 
from the surface observations and upper-air data. These 
parameters are defined as follows (all temperatures in degrees 
Gellstus, )nelghts in meters, wind speeds in knots and pressures 
Prem. liibars) : 

Ve ue: Pie wtemperatunre at the base of the ianverston ae 


OmeemnverSdon exiSts, Otherwise zero, 
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Bee ciee ea tutemate tie ton, Cr “the inversion if 
ane imversSlon exists, otherwise zero. 

Enews Strength of the inversion, measured as the 
difference between the temperature at the top 
of the inversion and the temperature at the 
base of the inversion (TT-TB=SI), if an inver- 
sion exists, otherwise zero. 

PVC dmemoreme base Of the anVersion tf an 
PMVwerctmonwexists, otherwise zero, 

Cie nelGiemoemsune top Of Ehe inversion a4 Jan 
inversion exists, otherwise zero. 

bnew theokness of the inversion, measured as the 
aqifference between the height of the top of the 
inversion and the inversion base (HT-HB=THK), 
if an inversion exists, otherwise zero. 

emes ct ecction from which the wind is billowing 

SE the 950-mb level. 

the wind speed at the 950-mb level. 

the temperature index, measured as the differ- 
ence between the temperature at the top of the 
Pigs mem Orerie SuUriace dry—-bulb> tempenature 
if there is no inversion) and the sea-surface 
pempenature (TE tor TDB)=SST=TI). 

the surface wet-bulb depression, measured as 
the difference between the surface dry-bulb and 


wet-bulb temperatures (TDB-TWB=WBD). 
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mie DPD: the dew—-point depression, measured as the dif- 
ference between the surface dry-bulb and dew- 
point temperatures (TDB-TDP=DPD). 

Hee WET: the surface moisture index, measured as the 
difference between the dew-point and the sea- 
surface temperatures (TDP-SST=WET). 

3 TDB: the surface dry-bulb temperature. 

14. TWB: the surface wet-bulb temperature. 

15. TDP: the surface dew-point temperature. 

16. SST: the sea-surface temperature. 

ey RH: the surface relative humidity. 

es a SLP: the sea-level pressure. 

foo FS: the fog strength; an integer number from 0 to 
4 indicating the number of three-hourly surface 
observations during a twelve-hour period (from 
Of OO=1600 PST or 1600-0400 PST) reporting fog. 

The notation used to express the parameters and their 
relative day and time within the three-day history period is 
reese tie wee Tne pre-subscript denotes the day during this 
Powuod a Wwikhle the post=subscript denotes the hour of observa-= 
meonmcOocOCsor 1600 PST), Examples of this notation are as 
Powwow si: 

ies KG = the height of the top of the inversion, taken 

Puonmenewmntrst.day Of history and the 1600 PST 


report. (This would be the most recent observa 


tiem an the history period.) 
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BF 04 = stile StuUenget, Of the inversion, taken from the 


Eien aucddyemens History andthe O4007PST report. 
(This would be the oldest observation in the 
Nisteuy serLoa.. ) 

In addition to the nineteen basic parameters described 
Pwo, 22—-, Anaeco-nour time differences in each of 
these variables were calculated using the three days of 
"history" associated with a fog occurrence. Because of the 
BMDO7M program Festrererroneot 80 variables there could be 
computed only 76 (of a possible 95) 12-hour differences, 76 
24-hour differences and 38 48-hour differences from the 
total list of 114 (1.e., 6x19) parameters for each day. 
These differences can be shown schematically as follows: 


12-hour differences 


| as (19 variables) 


Seaays of 0400 
history as- 3 1600 sk 
sociated aponee i 
with a fog 5 1600 ys 
occurrence 0400 " 

5x19=95 possible differences 

24-hour differences 
1600 (19 variables) 

11,0400 ‘ 
2 1600 - 


0400 si 


4 16 OC : 


jo Xe 


0400 a 


4x19=76 possible differences 
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48-hour differences 
1600 (19 variables) 
0400 " 
1} 1600 " 
2 0400 " 
Loo : 
0400 a 
2x19=38 possible differences 


Examples of the notation used to represent these differences 
eae : 


A = ther first (1600-0400) 12-hour difference 


eacao4 
in the temperature at the base of the in- 
version. 

A.TB = the second (1600-1600) 24-hour difference 
in the temperature at the base of the in- 
version. 

A TBO = the (0400-0400) 48-hour difference in the 
temperature at the base of the inversion. 

The time differences in the nineteen basic parameters 
are considered to be important since they allow for the in- 
fluence of persistence, advection and other physical changes 
iMimene Variables. Thus, the analyst can, at least, theoreti- 


Solvyspree yvect “trends” (time differences) of the parameters 


ahead in time for the purpose of developing a forecast. 


Cx DISCRIMINANT ANALYSIS~ 
The use of linear discriminant analysis in this study has 


been of fundamental importance in the successful development 


eTaken PeGmmcindmator e235 (written by William RK. Kleeka) 3or 
Statistical Package for the Social Sciences, Norman H. Nie, 
Pale ee Ghaw—-Hill, 1970. 
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GeeporassSification equations for marine fog prediction. The 
objective of the analysis is to statistically distinguish 
between two or more categories of cases. In this research 
these categories (see Table I) were determined by the occur- 
rence or non-occurrence of fog during three successive days, 
and whether fog was reported during the current 24-hour 
period. Thus, the groups “fog" and “no-fog" were established 
as previously described. 

In order to distinguish between the fog and no-fog 
groups, the nineteen initially selected parameters and their 
time differences were used as "discriminating variables" to 
Measure characteristics on which these groups are expected 
to differ. The mathematical objective of discriminant analy- 
Srsersommomweilgnt and linearly combine the discriminating 
variables in some manner so that the groups are forced to 
be as statistically distinct as possible. In other words, 
one hopes to be able to “discriminate" between the groups 
in the sense of being able to tell them apart. 

Discriminant analysis attempts to accomplish this by 
fmmine sone OL more linear combinations of the @iscriminat- 
imeoevarLabples. These "discriminant functions” are of the 


form: 


ib =e: +e .V. +c _.V. + 2.26 tec Vv (1) 


where De Pater alscraninmank Score for group mm, the ons 
qeewtliemecoert ficients with oc; being the “constant, enc meme 


Vsmerarcer the raw values of the n discriminating variables 
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used in the analysis. There is always a separate equation 
Roane dcmrOUuDweEOoperully, the discriminant scores (D's) 
bOimmthe cases Within a particular group will be fairly sim- 
ilar. Nevertheless, the functions are formed in such a way 
as to maximize the separation of the groups. Once the dis- 
criminant functions have been derived, the pL imary research 
Spyectave Of Classification can be pursued. 

Erassification 1S the process of identifying the group 
to which a fog/no-fog case is most likely to belong when the 
Only information known is the case's values on the discrim- 
inating variables. The classification procedure involves 
the use of a Separate linear combination of the discriminat- 
Migmevatiables for each group. These produce a probability 
of membership in the respective group, and the case is as- 
Signed to the group with the highest probability of occur-= 
rence. One problem with this type of classification proce- 
dure is that the rule of highest probability defines a very 
Searetecmyiding tines) A «51 probability of a fog occurrence 
versus a .49 probability of no fog would necessarily lead 
Mememrogwelassification, yet the situation is not really all 
thiatmaclLear . Mebavyesian adjustment Of the posterior proba— 
bilities is often desirable when the costs of misclassifica- 
PEMA tOmcerealin Groups are high, or when the groups are of 
Grosso me, sattererent sizes, In this study the posterior proba- 
bilities were not adjusted. 

During the analysis phase of this study the primary use 


GQewelacsccmrircation was to test the adequacy of the derived 
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Sisecrininant functions. By classifying the fog and no-fog 
cases used to derive the functions in the first place and 
comparing predicted group membership with actual group mem- 
berchip, One Can empirically measure the success in 
SiscElmination by observing the proportion of correct 
Slassiitcations. 

icra ho, eiprogmaneuced In this research 1s 4a Stepwise 
linear discriminant analysis routine which can accept a 
maximum of 80 discriminating variables. Since there are 
generally more variables available than are necessary, the 
stepwise analysis procedure is very helpful. It begins by 
choosing the single best-discriminating variable as deter- 
Mea by the Selection criterion. A second discriminating 
variable is selected as the variable best able to improve 
the value of the discrimination in combination with the 
first. The third and subsequent variables are chosen in a 
cam Lax Manner according to their ability to contribute to 
further discrimination. At each step, variables already 
selected may be removed if they are found to reduce discri- 
mination when combined with the more recently selected 
variables. This process continues until the remaining vari- 
ables are no longer able to contribute significantly to 
Pwenewr ad) Semi minatlon. 

The stepwise analysis procedure just described has sev- 
GHalmitrest Criteria which control the entry or removal of 
Wie cmriom the discrimination equation. Ultimately, 


the selection criterion is determined by the F ratio or 
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Mapeebue. the FE value ws the likelihood ratio of equality on 
Pitemtest VYalbtiable over all the groups (two groups in this 
Boas je Cl Ven Enon Ga strrbution produced by the variables al- 
ready entered. In other words, the "F-to-enter"™ value is a 
miter or thiematattoctical Significance of the amount of cen-— 
troid separation added by this variable above and beyond the 
separation produced by the previously entered variables. At 
each step in the BMDO7M program there is a Similar test made 
of all variables already selected which have an inclusion 
Pweimoc One (riven iehuston—-level default option Of one was 
used throughout this study). Here, the test is whether the 
Deattwecular Variable Still adds a significant amount to the 
group Separation, given the other variablesS now in the equa- 
tion. AS more variables are chosen, it 1s possible that 
some of those entered earlier will no longer be contributing. 
This occurs because the information that they contain about 
group differences is now available in some combination of 
the other included variables. The contribution is measured 
by the "F-to-remove", which is an F test of the discrimina- 
Prone power currently added by the variable in question. Tf 
a variable already entered has an F-to-remove value less than 
one (again, the default option was used), it will be removed 
PrOnernieweduation. 

The BMDO7M program computes unsStandardized coefficients 
Pm emetassttrcation LUnctions, thus making analysis and 
Piemonte tatton Of Such coefficients difficult. The unstand= 


ardized coefficients do not report the relative importance 
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of the variables since they have not been adjusted for the 
Measurement scales and variability in the original variables. 
However, the unstandardized discriminant function coeffi- 
events Can be Gasily wtilized to find the fog/no-fog proba- 
bapeetiy. In Order tO do so one first computes the 
Geseriminant score by multiplying the coefficients by the 

raw values of the associated variables, summing them together, 
and adding a constant to adjust for the grand means. This 
process 1S represented by Eq. (1). Then the posterior proba- 


bility of case k having come from group m is: 


exp (D ) 
+ = 
m,k Pe PID Ve ae Put EON sel 


where Ne aerecgmor no=-Log day (f or n£) 
oe thiemorEtor probability of group m 
D ea miemitscriminant Score Lor case k of group me 
f 


iMimenmmsectway stne prior probabilities of belonging to either 
the fog or no-fog group were taken to be equal, so P,, was 
Woe SOmMe aLSCYriminant analysis programs allow for automa- 
EPeemadgustments tO the pricer probabilities on the basis that 
SUcwmprebaptttEl1eGs are proportional to the number of cases 
in each group, i.e., cases will be more likely assigned to 
aplargqer group. The BMDO7M routine does not have an auto- 
Matic adjustment feature, but manual changes in the prior 
probabilities can be accomplished. 

A wealth of statistical information can be derived from 


the BMDO7M computer program--much more than this researcher 
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Pogwrrea se lie inzOrmatlion generally used in this study con- 
sisted of: 
1. Group means and standard deviations 
Cage At each step: 
a. Variables included and F value to remove 
Deva iables Wore included and F value to enter 
3. After each of the first seven steps and after the 
Maseestepmetie classification functions were printed 
4. tem posterior probability of coming from each group 
was computed for each case 
ae Summary table. For each step of the procedure the 
following were utilized: 
a. Variable entered or removed 
love F value to enter or remove 


ee Mumoer OL Variables included 


i ANALYSIS ROUTINE 

After the data had been organized by location, season 
and fog category, the job of analysis was finally initiated. 
Since considerable fog research has been done in the San 
Diego area (see Section II), this study began here. 

POrmeadcietlocatilon and for both the “dry" and “wet" sea= 
sons, the fog and no~fog cases of categories 1 and 2 (and 
later 3 and 4) were combined and read into the BMDO7M rou- 
tine. It is important to remember the program restriction 
SieaMaxrmumsOrecO GisScriminating variables, since this lam- 


itation influenced the analysis methodology. Associated with 
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eaommemog or NOo-fLog occurrence were 3 days (six reporting 
perwods) Of History, Listing a total of 114 variables (6 
reports x 19 parameters = 114 variables). Also recall that 
there were the 12-, 24-, and 48-hour time differences of 
these variables which were available, thus adding another 
POGeDossibile fog discriminators. However, for any one com- 
puter program no more than 80 of these variables could be 
Wie Zed . 
iieOnder tOetest@Ene discriminating potential of all the 
possible variables, seven separate analySisS programs were 
completed for each location, season and category pair (1 and 
See emacs Lhese, programs differed from one another by 
virtue of the variables being read into the BMDO7M routine 
and are described as follows: 
1. The 19 initially selected parameters for only the 
tHHietewordays Of history, but both 1600 and 0400 
PST reports were used. 
2% theme basic parameters for all 3 days of history, 
Timneeeoie lOO r ST repores were used, 
MaeuncCmlowodsloe DPatameters for all 3 days of history, 
Mitmenityvesene O2OQO PST neports were used. 
4. Only the 12-hour time differences of the 19 basic 
parameters were uSed (four 12-hour differences). 


Doe tie s2z4—-hour time diiferences for only the 1600 PST 





reports were used. 


Speeelne 24-hour time differences for only the O400 PST 





reports were used. 
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Zee wicLo-WOwlm tame ditferences for both reporting 
periods were used. 

Each of these computer programs executed the discriminant 
analysis routine with the available variables and printed 
a summary table listing the discriminating variables in 
Gtieder Of entrance and the classification functions. After 
Soret ing gttese "jobs tor both seasons (and for a given fog 
category pair), the best seven or eight discriminating 
variables were selected from each program and put ina list. 
Any of these variables which may have been duplicated be- 
tween, say, programs #1 and #3 (see list above), or between 
wet and dry seasons were discarded. Finally, the “best" 80 
discriminating variables were computed ina "composite" 
Maaogram Lor final analysis. This composite program was 
executed for both seasons and resulted in the final classi- 
Puecdetonetunetions £LOr a given location and fog category 
pair. Then the entire process wasS repeated for the other 
Meg category pair (3 and 4). ta Sopot Oar eae i locatron. a 
minimum of 32 discriminant analysis programs were executed 
MimeGmaaem tow build the ELEinal classification equations. When 
uSing these equations one can calculate the percentage pro- 
bability of fog/no fog over a 24-hour period beginning at 


1601 PST on the current calendar day. 
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V. ANALYSIS RESULTS 


A. INTRODUCTORY REMARKS 

The analysis results of this research are presented 
largely in tabular form. For each west coast station studied 
Pictemarc tChreenarnrenent types of tables; The first type 
Shows the number of usable fog and no-fog "cases" (sometimes 
called "foggers" or "no-foggers") for each analysis category. 
Ausable fog (no-fog) “case”" is simply a fog (no-fog) day 
for which there are a complete list of parameters. The second 
type of table presents the classification functions, the 
variables utilized, and their associated F values. Recall 
that the F value (F ratio) shows the relative importance of 
the variable in question to the existing equation. The F 
values shown in the tables are those for the final classifi- 
Sietonweguatione, Thus, in general, the larger the F value 
Ccetotiverco unity) £LOr a given variable, the greater that 
variable contributes to the discrimination (centroid separa- 
tion) between the fog and no-fog groups. The third type of 
table shows the season and er esory FOrRargLVGhw Loc at rem, 
BiGmeme peEtcentage of fog and no-fog cases correctly clas= 
Strteamey) the Giscriminant functions. 

It is important to remember that the stepwise dis- 
Cramtnont analysis procedure results in an optimal set of 
Variables being selected and that the initially-best dis- 


niet ing variables are chosen first. Therefore, in the 
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Pees PpLeEesenting the Classification functions the variables 
Se eotea ti Order of decreasing importance according to 
tame rat entry. Where a variable had been entered and was 
subsequently removed there appears a remark indicating the 


teepminm the program that the variable was deleted. 


B. SAN VOLE GOSRESULITS 

Repeal sang 1 Lia chow the number of fog and no-fog 
cases for the San Diego dry and wet seasons, respectively. 
During the dry season the fog cases comprise only 15.5% of 
the total number of cases, while during the wet season over 
28% were fog days. 

Table IV presents the fog and no-fog classification func- 
tions for the dry season, categories 1 and 2. The F values 
piarGate that the two most important discriminating variables 


in the final equations are the temperature index ( ) ®aned 


1704 
the inversion thickness (, THK) 4) taken from the most recent 
O400 PST observations. In this case the mean 1ito4 for the 
foggers was approximately 4.0 C, while the no-foggers had a 
Mesnonrmasoute).6 C. This implies that a relatively large 
temperature difference between the air temperature (at the 
Hmeacecmeare inversion top) and the SST prevailed just prior to 
most fog occurrences. This result generally supports Leipper's 
(1948) temperature index range values found in San Diego. 

The mean values of the inversion thickness for the fog and 


MOO gmGimsces aqittered only slightly (385.7 m versus 362.3 5m); 


DWE nemmarcKMessS was greater for the fog reports. 
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Mien olassiiication CLunctions for the wet season, catego-= 
ries l and 2, are given in Table V. Although the moisture 
index (, WET 


) and the inversion strength ( ) were en- 


ie D4 


tered as the first two variables, they were subsequently re- 
moved. When they were later re-entered they made only a 
BMtatiecOmer Out toOlmto Lhe overall discriminating power of 

the final equations. In this program the ultimate best three 


variables are the wet-bulb depression (, WBD );, thesheichtac 


16 


the top of the inversion ( HT yy) and the temperature index 


3 


( ) at 1600 PST. The mean wet-bulb depression for the 


eer 1 6 
no-foggers was almost 1.3 degrees greater than that of the 
Begerepores. 9 Lhe Ehnee-day-old HT value was, on the average, 
81 m higher for the no-fog cases. The TI mean values dif- 
Homo yeanoOomte Lol € (35.3 versus 2.2), and again this value 
was usually larger for the fog reports. 

Table VI shows the classification functions for San Diego 
during the dry season for categories 3 and 4. Recall that 
Gategories 3 and 4 differ from 1 and 2 in that they have a 
Mistory Of fog. Furthermore, categories 3 and 4 should give 
some clue as to the parameters of interest in fog dispersion. 


) 


The F values once more identify the temperature index (,TI,, 
as being very important to the final equations, and surpris- 
Pig, eco 1s the three-day-old wind direction (,DIR,,) > As 
was seen before, the mean TI was significantly higher for 
the foggers than the no-foggers.j. In this instance the dif- 
PomomcecmuasmabOUtL 2.5 degrees (i1.e€., 5.5 C for fog versus 


BUme mE OLeNO-fOg). in San Diego the temperature index ap- 


pears te be an important fog/no-fog discrimination parameter. 
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The 950-mb wind direction mean values indicate that 
three days before a fog occurrence the wind was southeasterly 
(about 126°), whereas prior to no-fog days the wind was gen- 
Ssaaulyet COMstne ~—SOUuthwest (207°). It should be noted, how- 
ever, that the standard deviation for these wind directions 
was about 120 degrees. Therefore, solid conclusions as to 
mie sSsigniticance OfmeEniSs result cannot be made without further 
muriysts On a broader Gata base. Other important discriminat— 
ing parameters shown in Table VI are a 48-hour change in the 


1600 PST wet-bulb temperature (ATWB ), a 48-hour change in 


16 


the 0400 PST wet-bulb depression (AWBD),), and a 12-hour 


yc 


change in the temperature index (ATI 
J P a omens 
The San Diego wet season classification functions for 
categories 3 and 4 are presented in Table VII. The most 
capable discriminators in these equations are the wet-bulb 


), the inversion-base height (.,H ), and 


depression (, WBD l Bis 


IEG 
the fog strength (_FS lemeali measured at 1000 PST: The 


eevee hs 
Statistics for the wet-bulb depression mean values are nearly 
identical with those mentioned for categories 1 and 2 during 
the wet season. The no-fog reports averaged over 1.3 degrees 
hagher than those of the fog cases (4.3 C versus 3.0 C res- 
PeG@eveia) . thas result appears to be consistent for all 
categories in the San Diego winter season. The mean values 
Sreeiewtnverston—-base height exhibit dramatic differences be- 
tween the fog and no-fog cases. in Log Situations the’ bace 


heights were about 168 m, whereas the heights for no-fog 


SPueuatiemo were 536 m. Thus, as the fog regime begins to 


47 





“rcciparte, the base of the inversion rises. This is in 
agreement with Leipper's model for San Diego winter fogs 
(Leipper, 1948). 

ive wteOguSotrengtn (FS) parameter can be a significant 
meserciminator Only in Categories 3 and 4 where there is a 
[mist Oory Of Loge In this case, during the wet season the 
mean FS value for 1600 PST was nearly 0.88 for fog cases 
Amcmeer zo LOL NO—-fOog reports. This simply suggests that fog 
meesistence is Likely to occur in winter. That is, fog is 
ice ly tO occur in one of the three-hourly surface observa- 
tions in successive twelve-hour periods. 

Table VIII shows the San Diego fog and no-fog discrim- 
ination percentages for the seasons and category pairs. The 
Sscriminating Capability of the category 3 and 4 equations 
is better than those for 1 and 2, especially for fog classi- 
fication. Overall, the weighted average fog discrimination 
percentage was about 82.2, while the no-fog percentage was 


mearly 81.3. 


OR OAKLAND RESULTS 

The occurrence of fog in Oakland's summer (dry) season 
is indeed infrequent. For the five-year data period of this 
study there were only a total of 32 fog days during this 
season, as shown in Table IX. This amounts to less than 
POUm@m~Deomcemt wear tChe total cases. Fog occurrences in the wet 
season (shown in Table X) were more common--about 24% of the 


Gases reported fog. 
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Boole exXispresents the fog and no=-fog classification func-— 
mons £Or the Oakland dry season, categories 1 and 2. As the 
F values reveal, the best-discriminating variables in the 
momab equations “are the 24-hour change in the 1600 PST in- 


memsion thickness (A, THK je ra) SimMillJare time ditterenee niet ne 


IG; 


relative humidity (A, RH Yeeand the second (0400-1600) 324— 


16 
hour change in the temperature index (AjTIQ,)- It is in- 
m6esting EG note that the most important variables in this 
equation are all 24-hour time differences. With the exception 


of one wind speed variable (SPD ), the nineteen basic para- 


Ins 
meters themselves contribute very little to the fog/no-fog 
Greserimination. 

The 24-hour thickness change mean values showed that 
there usually was a decrease in the 1600 PST inversion thick- 


ness just prior to a fog occurrence. The A, THK average 


16 
for the fog reports was about -132 m, whereas the no-fog 
value was 2.6 m (a negligible increase). Although the stan- 
dard deviations for the two groups average about 265 m, the 
means clearly indicate a reduction of the inversion thickness 
which was an unexpected result. Since this season contained 
Only 25 fog reports, analysis of a larger data set seems 
necessary in order to be certain of this thickness relation- 
Slip. 

Aaoenmer Surprising result in this analysis was that the 


2a=nowre relative humidity change (A, RH ) for the fog reports 


16 
waS an average -0.18%. Although the mean change is very 


Stewiewone would normally expect a moderate increase in the 
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Po oteve nuMmrdity prior to fog. A for the no-fog 


i 216 
cases was almost zero, and the standard deviation for both 
Geomps exceeded 2.0. 

The mean values of the 24-hour change in the tempera- 
ture index (A,TIOy) reflected an increase (almost 1.5 C) 
for the fog cases and a very slight decrease (-0.06 C) for 
the no-fog cases. This result appears to be consistent with 
San Diego observations; several days prior to fog develop- 
ment the inversion intensifies which generally leads to an 
increase in the temperature index. 

The classittecation Functions for Oakland's wet season, 
categories 1 and 2, are given in Table XII. Here the filmuse 
two variables entered in the stepwise discriminant analysis 
Memarnea the best variables in the final equations. Both 
of these parameters were from the most recent 1600 PST ob- 
servation and are the temperature at the top of the inver- 
), and the relative humidity ( 


Son. «( RH, ¢) > The means 


ines 1 


eeeche TT values were 8.7 C for foggers and 5.5 C for the 
no-foggers. From this result one might infer that the in- 
version is either thicker or more intense, or possibly both. 
The RH mean value for the fog caseS was almost ten percent 
higher than for the no-fog cases (73.5% versus 63.63%). 

Of tertiary importance as a discriminating variable is 
the first (0400-1600) 12-hour change in the temperature of 
ENeCwInVvVersion base (4, TBO 4-16): The mean values for the fog 


Bawa —hOGmeases were approximately 3.2 C and 1.4 C, respec=- 


tively. Although the inversion base showed a 12-hour warming 
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tendency for both groups, the average increase was more than 
twice as great for the foggers. The reason for this warming 
Eemaency at the base of the inversion is not fully under= 
SLood by this researcher, especially since this was over a 
mewour perrod (G400=-1600) during which diurnal cooling us-= 
Maely OCCUurs. 

Vabbhe XIII shews @ebe classification functions for Oak- 
land Guring the dry season for categories 3 and 4. As the 
Pevaluecs @magtcate, the three most capable discriminators in 
Mlese equations are the fog strength (.FS._), the 950-mb 


2 16 


wind speed (SPD ), and the wet-bulb temperature (, TWB ee 


16 16 


Curiously, all three of these variables were taken from the 
MoeOerPsST reports which were two days old. It is important 
to note that there were only seven fog reports for this sea- 
son and category, compared to over eleven times as many (78) 
no-fog cases. 

heewestangqly, the fog strength mean for the fog cases 
(0.14) was less than half of the mean for the no~fog cases 
reo) eee oes gust the opposite that one might expect 1£ 
fog persistence had occurred. Because the fog and no-fog 
cases were so disproportionate, perhaps some adjustment in 
the prior probabilities reflecting this imbalance would 
moderate or change this result and remove a potential sta- 
Prot Lcals bias. 

The wind speed parameter showed less than one knot dif- 
ference between the mean values of each group--3.9 knots for 


boggqens sand 4.8 knots for no-foggers. The wet—bulb 
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PeMPeianelhe averages were about 16.0 C and 14.3 C for the fog 
and no-fog cases, respectively. Once again, a larger data 
sample of reports is required to prove whether these differ- 
ences are really significant. 

The Oakland wet season classification functions for 
categories 3 and 4 are presented in Table XIV. The best dis- 
criminator in these equations is the temperature at the inver- 


Sion top, measured at 1600 PST (,TT )--the same as for the 


16 
category 1 and 2 analysis during this season. The mean 
values for this parameter indicate nearly a 4.3 degree 
Spread between the groups--about 10.0 C for the foggers and 
5.7 C for the no-foggers. The identical inferences may be 
made since the tendency was equal to that of the category 1 
euic) 2 result. 

TvOmOotncueLmportant discriminators in this program are 


the 48-hour change in the 1600 PST 950-mb wind speed (ASPD.__), 


WS 


amas the O400 PST height of the inversion top (HT Ja ihe 


04 
average wind speed change was only a one-third knot decrease 
for the fog cases and a 0.6 knot increase for the no-fog 
reports. Nevertheless, in this data set the nearly one knot 
difference was found to be statistically significant by the 
EMbOyMaeprogram., The mean values for the height of the tin- 
version top for foggers and no-foggers were 566 m and 637 m, 
Cespective ly .« Although the difference between the group 


means was only 71m, the tendency for a lower near-surface 


iver stronmpitor to fog occurrences is evident. 
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Mawiwe “Aye shows Ene Oakland fog and no-fog discrimination 
Pemeemreages Lor the Seasons and category pairs. The high 
Peceaiminatton Seore for the dry season, categories 3 and 4, 
Maoeeaare ly sauce COSthe unusually small number of fog and no= 
fog cases. If the BMDO7M program had been allowed to con- 
mime Unter! the F=to-enter value dropped below 1.0 (as was 
done with all other equations), the resultant equation would 
have had 30 variables (after 48 program steps!) and could 
memieve a vpierfect” discrimination score. However, thirty 
variables was considered to be an excessive number for this 
small data set and the program was limited to twenty-one 
Steps which included the final fifteen variables shown in 
feoremrriics The Overall weighted fog discrimination percent=— 


age was above 77.4, while the no-fog percentage was 76.8. 


D. CULTEBDAYULE RESULTS 

Whereas the occurrence of fog in Oakland was a rarity in 
summer, Quillayute, Washington has fog as a fairly frequent 
Wetter tnhhougiout the entire year. In fact, for the period 
@ieschis research exactly 50% of the dry-season reports and 
nearly 54% of the wet-season reports were fog cases. The 
number of fog and no-fog caSes for these two seasons are 
shown in Tables XVI and XVII. It is likely that in this area 
a significant percentage of fog occurrences were associated 
Vilg@li@merOmbal activity, Since most mid-latitude storms track 


across the Washington coast. 
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Pacey ttt epresents the fog and no=fog classification 
PGmetlons tor the OQuillayute dry season, categories 1 and 2, 
There are many variables with relatively large F values, but 
the three best discriminators in these equations are the 
iescerecent 24-hour change in the 0400 PST temperature index 


(A, TI ), the second 12-hour change (0400-1600 PST) in the 


04 


temperature index (ATI i, ana the 950-mb wind direct uon 


04-16 


measured at 0O400 PST (, DIR ae Also notice that the “furccte 


04 
two variables entered in the stepwise discriminant analysis 
(A, THK), and AWET, ¢) were subsequently removed. 

The importance of the temperature index parameter and its 
time differences iS again demonstrated by the fact that they 
are the top two discriminators in these equations. The mean 


values for A,TI Were aqoout 1.4 ¢C for Loggers and 0.6 7Cerer 


04 
no-foggers. Thus, the 24-hour change in the temperature 
bmieiex was Significantly higher for fog situations. This 
result is likely to be caused by increasing temperatures at 
the top of the inversion which is associated with a strength- 
ened inversion. The mean values for ATI 54-16 were about 
moe em ance 3.6 C Lor the fog and no-fog cases, respectively. 
Maemeugh the BMDO7M program found some statistical signifi- 
cance in the group means' difference, the standard deviations 
were fairly large (about 2.4 for both groups) and very few 
conclusions can be made about this parameter. 

The 950-mb wind direction averages were about 214° for 


Peoqmeewoarts and 186° for no-fog cases. Again the standard 


deviations were large (nearly 115° for both groups), but the 





lemming Girection prior to most fog occurrences had a com- 
ponent from the west. This implies that the fog was generally 
advected landward by the prevailing southwesterly flow. 
ftewelassitwedemon Functions for Quillayute's wet season; 
eauegories 1 and 2, are given in Table XIX. The first var- 


iable entered, the 1600 PST 950-mb wind speed (, SPD ), was 


16 
the best-discriminating parameter in these equations. The 
mean wind speed for the fog cases was over 9.0 knots, while 
mme no-Log mean was 5.7 knots. The cause for the signifi- 
cantly higher wind speed associated with fog occurrences is 
not clear. The supposition that these speeds are related to 
pre-frontal fog appears to be disproved by two other impor- 
Panerai scriminators. 


The two-day-old wind direction taken at 0400 PST (DIR ) 


04 
was the second best-discriminating variable. Surprisingly, 
mie mean wind directions for both fog and no-fog groups were 
from the southeast gquadrant--the values were approximately 
[Meo and 125°, respectively. The third important variable 


was the sea-level pressure measured at 1600 PST (,SLP a 


16 
The group means showed a Significant pressure difference of 
aeouer 2.7 imb. The average SLP for foggers was almost 1019.6 
mb, while the no-foggers had an average of 1016.9 mb. This 
researcher was unable to find an adequate explanation for 
this last result without supplementary information. 
Wimieomeeweschows the classification functions for Quatlla-— 


yYyuce Guring the dry season for categories 3 and 4. The F 


Voloeouwindiwcate that the first-entered variable, the 1600 





PST relative humidity (,RHig) is by far the most important 
aeseriminator in these equations. The mean relative humidi- 
eiecrOr the fog and no-fog cases were 73.5% and 60.4%, res= 
pectively. it appears logical that as fog approaches the 
Mette romeene relative humidity will increase. Tie seks 
Spread between the average relative humidities provided easy 
discrimination between the two groups. 

Mecemetie: —vud~iaples im lable XX stand out as signifi-= 
cant in the classification functions. These are the tempera- 
ture index /( 


), the 950-mb wind direction (, DIR )}, and 


Te ite 


the moisture index (, WET ), all derived from the most recent 


iG 
1600 PST observations. The TI mean values were -2.8 C for 
meqqers and about -Il.7 C€ for no-foggers. The fact that the 
fog cases had a larger negative TI is rather surprising since 
one would expect that the closer the air temperature is to 
memos, tme Greater the chance of fog. The negative TI in-=- 
dicates that the air temperature is usually lower than the 
Sol. 

The mean wind directions for the fog and no-fog groups 
PomemoOoth trom the Southwest (237° for foggers and 231° for 
no-foggers), and showed Such a minor difference that one 
wonders why it was chosen as a meaningful discriminator. The 


moisture parameter (, WET ) seems to have greater theoretical 


16 
meaning in this program. The mean moisture values were 
oe ronmeately —=—Seo-e for foq cases and -9.9 C for no=-fog 


reports. Recall that the moisture index is the difference 


between the dew-point and sea-surface temperatures. In the 
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QOuilliayute dry season the dew-point temperature is nearly 
always well below the SST. However, during fog situations 
the difference is smaller, which implies a greater likeli- 
iIMgod "Or Log Lormation. 

Timresourgitlayute wet season classification functions for 
categories 3 and 4 are presented in Table XXI. In these 


equations the relative humidity parameter ( ) again 


ine 
Seands out as the most capable discriminator. As before, 
the difference between the groups! averages is about 13%, 
but in the wet season the relative humidities are higher-- 
Slmestlore0s Lor Eoggers and /4.2% for no-foggers. 

The second- and third-best discriminating variables are 


the 1600 PST temperature index ( ) and the 0400 PST wet- 


1 16 
bulb depression (,WBD,,)- ine “Ti stetrstics showed farrily 
large negative mean valueS--a -7.7 C for fog cases and a 
-8.1 C for no-fog reports. During the winter season one 
Biewlommexpece Larger negative values for this parameter in 
the colder, more northerly latitudes. Nevertheless, the 
Slightly smaller TI for the foggers indicates a greater 
probability of a fog occurrence. The wet-bulb depression 
means values showed only a minor spread of 0.1 degree--the 
foggers averaged 0.54 C while the no-foggers were about 

0.64 C. Although the difference was small, the discriminant 
analysis routine found that it made a significant contribu- 
PUGnMmnimtnese clasSifiication functions. 


taste wil presents the Quillayute fog and no-fog dis= 


crimination percentages by season and category pairs. The 
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Memiweeo tor the category | and 2 analyses are considerably 
Mecwcr tian Lor 3 and 4. This was probably due to the much 
pliliaiteresmumber Of cases in the 1 and 2 category. The over- 
Beiwergntead GdiISCrimination percentage was above 74.9 for 
muomtou Cases ana nearly 75.5 for the no-fog reports. 
iiemeuluiatave average Giscrimination percentages for 


Beeemtece Stations analyzed in this study resulted in 78.23 





memralilytog cases and ever 77.8% for all no-fog cases. 


Ee SOUE- LES Lp RESULTS 

PimeedeomaeomGest Ene GiScriminating capability of some 
Steche Classification equations three years of additional 
data were utilized. The Naval Weather Service Detachment, 
Asheville, North Carolina again provided the airways surface 
observations and radiosonde soundings, while more SST infor- 
Mation was extracted from the Fishing Information Bulletin 
Supplements. The test data covered the period from 1 July 
1965 through 30 June 1968--i.e., the three years preceding 
elie initial analysis period. 

Peindmiiysdaue to time considerations, only the San Diego 
classification equations were tested. The data were prepared 
by season and categories in the manner previously described 
(see Section IV, subsection A). The number of fog and no-fog 
Case omerst cas sy Category for the test data dry and wet Sea= 
sons are shown in Tables XXIII and XXIV, respectively. A 
testing computer program was developed which would form the 
pogqmancdane—-fog probabilities (see Eqs. (1) and (2)) and 


classify each case in the same manner as the BMDO7M routine. 
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iomeom en Ve Shows ethe percentage of correct discriminations 
Dmbetecesoulted from testing the San Deigo classification equa- 
Mincmenmm tine Independent data set. Although the discrimina- 
mEOMmmpoereentages were considerably reduced from the dependent 
analysis sample, the outcome was generally acceptable except 
Poem WewabyesoecasOn Log prediction capability. The reason £Eor 
mies tavlure 1S not known. The overall weighted fog discrim-— 
ination percentage for the San Diego test data was above 
68.4, Re the no-fog percentage was nearly 74.3. 

By comparing Ehe numbers of fog and no-fog cases in 
maabes Lit and XXIV; one will note that the relative propor- 
tions of cases are similar. Thus, it seems unlikely that un- 
Dewalee Group sizes could be a factor in the dry season's poox 
Geserimination capability. Perhaps better results could be 
obtained if the BMDO7M program prior probabilities were made 
Pmopertlonalycto the number of cases in each group. This 
Wonlla “inerease the likelihood of being assigned to the larger 
Groupe Another alternative would be to make a Bayesian ad= 
fiscemnont to the probabilities of group membership, especially 
when the groups are of grossly different sizes as was found 
Mroeotite scam Daego and Oakland. 

Paoimnable xx itl one can £ind that during San Diego *s 
pmeyescecasonmoenlty 11% of the total number of cases (for the 
EeStmecicinadapertod) were Log reports. During the wet season 
(com iasa lew xXxXi1V) about twice as many (22.5%) fog days oc= 
euured. Eeome these Statasties a forecaster might logically 


conclude that he would be correct almost 90% of the time in 


ao 





mmemanry season by never predicting the occurrence of fog! 
Indeed, persistence and climatology forecast methods may 


Mermatiticult to Surpass, but other prediction tools are 


essential. 
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Vale PROCEDURE FOR APPLICATION 
iS Sedat ts 
Shnocolr IT CATION FUNCTIONS 

Siew Otecthie Primary objectives of this study was to de-= 
fom we lassirication equations which could be utilized to 
Pec cast ete dally probability of occurrence Of marine fog 
Memelne three U. S. west coast stations analyzed. The clas= 
Perec ation functions presented herein provide a means for 
Meeoeasting the probability of fog/no fog for a 24-hour 
Pewwod beginning at 1601 PST. In order to calculate these 
probabilities the forecaster must have: 

a The station's 3-hourly surface observations and the 
upper-air soundings for the last three days. 

2. A reliable bi-weekly mean coastal SST near the sta- 
fe k@ Ti. (This is required to form the TI and WET parameters.) 

Seeeereeaesk calculator or a small programmable computer. 

The following step-by-step procedure 1S recommended for 
Mie orouer application of the classification functions. 

ieeNiote whetEher fog has occurred on any surface report 
in the last three days, Starting with the most recent 1600 
Bed Observation. bE enene has been suchwa, Log cecummence 
use the table for the current season (wet or dry), cate- 
Sories S and’ 4. PE Sener eeis no tistomy ot shod ye Gs c weme 
table labeled categories 1 and 2, for the current season. 

2. Using the surface observations, the upper-air sound- 


Micgcwemancdethne SST data, Find the numeric value of each of 
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tive Variables listed in the table selected in step 1. (See 
mre derinitions of the variables in Section IV, See econ 
Ey.) 

Seer bamwteme sata OL a desk calculator or a computer mul= 
tiply the raw values of the variables by their associated 
mag coefficients and add them together along with the con- 
Siant . Do the same with the no-fog coefficients. These two 
ame are the discriminant scores, D 


and Do (see “hq. SC) es 


f - 


4. TO avoid problems in exponentiating discriminant 
peeGtes which might be outside the allowable eae one may 
aqoerwtne Lollowing: 

a. Select the minimum discriminant score between 
bee and D 


de okays 


Dobe ract Lhe Minimum discriminant score from Dee 
ieemene result be called "DIFF". 
(1) ie Dinh ewisS Greater than or equal to 
ie 4.673, exponentiate DIFF. 
(2) Otherwise, set DIFF equal to zero and then 
exponentiate. 
Call the above result "PF". 
eumeeepeat step b using the no-~fog discriminant Score; 


D Mimomeces of De. Call this result "PNF”. 


ne fe 


EOrenneine sDOStEerLtOr probability of a fog occurrence 


by dividing PF by the sum of PF and PNF (See Eq. (2)). Then 


erie constant 174.673 is the largest power of e that 


Cicer oO COMpULeEr Can represent. if Using anotner seen 
Piltber Ehis constant will vary. 
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Gmempreobeability of no fog occurring iS unity minus the pro- 
Paiodlity of fog. 
Appendix A presents a generalized FORTRAN program of 


the above procedure. 
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Vii CONCLUSIONS AND FINAL REMARKS 


Dic phimMakty Goal Of developing linear classification 


equations which may be used for forecasting the probability 


of marine fog at San Diego, Oakland, or Quillayute, has been 


accomplished. The ultimate capability of these equations 


has not been fully tested. Nevertheless, it is hoped that 


these discriminant functions may serve some useful purpose-- 


if not for operational employment, 


mow Lurctcher study .« 


then perhaps as a basis 


One of the important by-products of the classification 


equations was insight as to which variables are the best 


nog, NO- tog GLsScriminators. 


Sreche Most capable discriminating variables, 


The following is a brief summary 


bas ted ach 


order of decreasing importance by location, Season, and 
Ga Leqgony Pair: 
SiN DLEGO 
DR veal et 2): 04 7 78S oa A,DIR, 6 _o4 AWBD) , 
Wee Ct); 1WBD, 6 «3HTy, itis ATB yy 
DRY 3-4): 17+ 94 30IR,, ATWB, | AWBD,, 
Vel wea) B51; peeis 1 16 ATI, « 
OAKLAND 
DRY (1+2): A,THK,, 4,RH,; ATI), ASST a_16 
PEE Eee + 2): tt 16 ee 4 TB Oy 16 3PIR,¢ 
De ero «SPD, 6 (TNE AN aya 
WET (344): ,7T,, ASPD,, ,4Tp, NS ISI ee 
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SULELAY UTE 


DRY (1+2): A,TIL, Py ae 1PIRG, ARH) ¢ 
WET (1+2): jSPD,, ,4DIRL, >SLP,, 45TT 6-04 
DRY (3+4): RH); i tie 1PIRi¢ WET 16 
WET (3+4):  ,RH,; wrtice “BDO, 4,SPD ¢_og 


Several conclusions can be made about the discriminating 
parameters shown in this analysis summary: 

i The PemMimeadeaurenanaex (TI) and several Gl a2tsetime 
differences are powerful discriminators, especially in San 
Diego and Quillayute. Since this index uses the SST aS an 
Pgehetepornt  , reliable coastal SST data are very important. 
ios May provide a clue as to why the TI is not particularly 
Seu tticant at Oakland. The SST data for this station was 
Sxtracted from SST charts at a point just west of the San 
Francisco Bay entrance, 


supposedly where the marine fog 


Might form and be advected shoreward. In ©etrospece.] per 
haps a bay water temperature would have been a better refer- 
emee fLor the TIL at this station. 

2. The most recent relative humidity measurement (, RH, 6) 


and its 24-hour change (A, RH ) 


16) 2Ppear to be excellent dis- 


Carirenating variables in Quillayute for Categories 3 and 4, 
and in Oakland for categories 1 and 2. 
the £600 5PST temperature 


SA In Oakland's wet Season, 


Beet emnvyerstion top (_TT..) is a useful fog/no-foeg pre- 


if 16 
Geetor < 


4. The 1600 PST wet-bulb depression is an important dis- 


Criminator in San Diego during the wet season. 
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ae Two discriminating variables which were hoped to 
be Of importance in the final classification equations were 
the strength of the inversion (SI) and the moisture index 
(WET). One will notice, however, that these parameters 
Beneraliy made only a relatively minor contribution to the 
feccriminating Capability of most equations. Several times 
the moisture index was entered into the regression equation 
Mmeesty Inaicating a large initial discrimination power, only 
to loose its importance as other variables were entered. 
Nevertheless, this researcher feels that the moisture index 
is still a worthwhile parameter and should be investigated 
iar urther research. 

The discriminant analysis approach is certainly not a 
new research technique. In the past, multiple linear regres- 
sion has been employed in numerous studies which require 
analysis of multivariate observations whose predictand was 
mMoem=-numerical. In the last decade, through the advance of 
Computer technology, tremendous progress has been made in 
mmerapid analysis of numerical and statistical information. 
Without the use of a high-speed computer and a sophisticated 
Prteatmnmnant analysis program the completion of this re=- 
Peamweimoy ed Single individual within a reasonable time frame 
Womladebenampossible. However, through the application of 
Pisce mpewerral computational resources meaningful relation- 
Stecwean be determined and an enhanced fog-forecasting 


capability may be evolved. 
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NEG EIE RECOMMENDATIONS FOR FURTHER STUDY 


The analysis results revealed several interesting and 
unexpected relations between various meteorological para- 
Merersesand the Occurrence Or noOn-occurrence of marine fog 
at selected locations on the west coast of the United States. 
imeeorcen to evaluate the full significance of these para- 
meters and attempt to improve the classification functions, 
the following recommendations are offered for future study. 

ime ebaccgetudata, base 1S needed inorder to obtain a 
greater number of fog and no~fog cases for discriminant 
Sar ysis. A Gistinct paucity of fog cases was noted in 
both San Diego and Oakland dry seasons. In such instances 
the BMDO7M program may generate classification equations 
which have used a few anomalous occurrences (which, however, 
Meee pe a significant portion of a small data set) for dis- 
crimination purposes. Such equations might not indicate 
the discriminating variables which are truly important to 
fogeforecasting. Consequently, a ten-year data period is 
considered to be the minimum for best results using the dis- 
CLaiminane analysis technique. 

Peace adjustments in the prior probabilities such 
Pieteeincce probabilities are proportional to the number of 
cases in each group. Further experimentation with variable 
aera obdiovlatiles 1s necessary to test their effect on 
Ptemsoteectluon Of the fog/no-fog discrimination parameters 


MmcCmencomeeiinat classification equations. 
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Emme coana, the List of discriminating variables to in- 
ede sviSibility, low cloud information and persistence 
parameters. Other wind levels should also be studied. 

Wane himMement ~Orere-definition of “fog” “to include 
only those cases where visibilities are reduced to less 
eeon chree mides (IFR conditions) may result in forecasts 
of greater operational importance. Some measure of rela- 
trve Log intensity, Such as "heavy" versus “light" fog may 
paso Cnhance the forecast. However, a very large data 
base will be required for such refinements. 

5. The use of better sea-surface temperature informa- 
tion is strongly recommended. Bi-weekly average SST data 
extracted from smoothed isotherm analyses are likely to be 
inaccurate. The SST data are important because a one- or 
two-degree temperature difference between the marine air 
and the underlying sea may be crucial in determining 
witether £Oog iS expected to form, or in measuring its inten- 
Setvee Reliable, daily SST anformation might make the tem-= 
perature and moisture indices (TI and WET) better fog/ 
no-tog Giscriminators. 

6. With the use of a larger data set, one should try 
analysis on four seasons of the year iImsteada Or Just twos 
Such a seasonal refinement may provide an enhanced fog pre- 
diction capability through the use of more representative 
Guasstuieadtilon Lunctions. 

foeserehop a method of stratifying the 24-hour fore-— 


cast into two 12-hour periods (or smaller intervals) in an 


68 





Siio@entO Indicate more precisely the expected time of fog 
meccurrence. 

8. Since Oakland is somewhat sheltered from the marine- 
Mim@erie uence by the coastal hills of San Francisco, data 
from a station closer to the coast should be used. The 
Surface observations from such a station could be coupled 
with the Oakland upper-air sounding data. 

Application of these refined analysis techniques to an 
enlarged data base would, hopefully, result ina significant 
improvement in marine-fog forecasting at the west coast 


Seatitons of interest. 
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CATEGORY HES TOR (o. DAYS ) CURRENT 24-HR PERIOD 





TABLE I. FOG AND NO-FOG CATEGORIES 
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CATEGORY ierpiigwe CLASSIFICATION 


50 FOGGERS 





531 NO-FOGGERS 






60 FOGGERS 






180 NO-FOGGERS 


TABLE I1. SAN DIEGO DRY SEASON FOG AND NO-FOG CASES 
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Gat iE GORy NUMBER & CLASSIFICATION 


69 FOGGERS 


407 NO-FOGGERS 


108 FOGGERS 
220 NO-FOGGERS 





ieee). tl. SAN DIEGO WET SEASON FOG AND NO-FOG CASES 
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GOCE aCclents ot the: 
FOG EUNCTION WO=FOG FUNCTION 








VARIABLE 

















FINAL F VALUE 


2241 16S 













ATWB -0.69860 -0.94497 










































































































































































































Io, -0.73689 =-1.07313 2376201 
A,DIR,¢_o, Or O05 0.00646 9.2365 
1THKy, O07 0012.5 0.00313 1 On? Si 
AWBD) , —Pe2332) 0.49745 S93 cic 
ADPD), OR 24705 =O. 27 8 5a, L529 
ASLP, | -0.36144 ~ O77 20 8.5494 
A,TWBOy OQ MiPo0.3 -0.09174 i 2007 
SIT a, 0200153 0.00081 BS 318 
ac G0 1062 Os Is 2 3.7899 
AU TDP, 6-4 1.19864 1.00863 292 Gal 
A,TI,, =O. 41794 -0.5445/7 1.8798 
ae =On2 4 5 5 -0.16593 4.2166 
PIR, , -0.00005 OS00234 Oh 6 (8 AISVe. 
,PIR,, 0.14416 Ovts es Da 2307 
Meo. 235922 6 2.24390 es SEO 
Constant -37.83261 = 3 5.4 ss 





TABLE IV. SAN DIEGO DRY SEASON FOG AND NO-FOG CLASSIFICA= 
POWeEUNCTIONS FOR CATEGORIES 1 AND 2 


13 





VARIABLE 


Coetratcients Of “thie: 
FOG FUNCTION NO-FOG FUNCTION FINAL F VALUE 


(Removed at step 16, & later re-entered) 


(Removed at step 20, & later re-entered) 
122204 O27 99060 239369 
On O08 7 > O 00260 4.6784 

-1.64342 =—le276o3 1 et 04s 
(Removed at step 22) 
Ze AO TSG 24 OO S012 3.4998 
=~1.Us4e7 =I 29956 A207 56 
Of 02139 5 OF 02075 3.4434 
0.00008 0.00134 LO. S558 
-0.07748 =0 boris T8624 
O26:07, 7-1 0.49668 4.1412 
-0.47578 -0.64106 Ax O329 
(Removed at step 25) 
Ze P14 Dae eS 515), S: CMP a5) (es: 
i O04 Le 27 535 1.9924 
OZOLe 75 On 0 Pes 2.4990 
0.00054 0.00094 3267 5 
(Removed at step 31) 
Of 93 22 O216020 dies) 0 his, 
= 2 A SNe SiS) -1.90545 2 leo 
P0573] O29 700 20290 
=o OOO 2 -2.83606 3.0366 
0 06976 O2. 19479 1. S140 
= 85104 = 1S ey ceals 2 Sel Jas 
Ieee sy syle. lees leo 265 
Constant -40.16452 -40.49783 





see TEGO WE Tesh ASOly POG AND NOP OG CLASse len 
MOM EUNCTIONS FOR CATEGORIES 1 AND: 72 


diol Ve 
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CoecEficients of the: 































VARIABLE FOG FUNCTION NO-FOG FUNCTION FINAL F VALUE 
FP 16 27: 5o6y3 126462245 3.4129 
3PIR,, 02472659 0343420 19.8636 
WET I ¢ (Removed at step 21) 
1itoa 2307572 -24.84428 222792 
F504 =—Je7 5544 — Oizo J 2230 
5TTo, 26.61739 265 70.6 G2 730 
AWBD) , a2. 25061 30.04062 10,2524 
ARH), 10795739 10.76887 Aves 374 
A,SPDi 6_o4 -38.42950 =360.25275 20 7 4a 
A.SPD O46 Da > 37 9 24232675 3.0246 
W 
AT Big (Removed at step 19) 
AoHB i ¢_o4 = (Ora 7e3;3 |! =O0:7,004 01 Vo 2383 
3°1i¢ (Removed at step 25) 
5SPD), J5.30409 75355360 4.1700 
prHKo, 0.41899 O42 0G / 3 43 52 
pee? og 221eeO2e30 Z2) 220645 324470) 
A, DIR), =—,12342 -1.12795 143105 
ATWB, . =o) oro 7 1 -97.94383 L22.096 
18D, ; 89.6240) 40.02596 6.3403 
4,7 416 6e.2396!1 68.77629 LOS 1973 
A,TDB,; 59.41071 59769287 120224 
3Tlig iO O05 3 7.34978 4.3841 
pePPo, -17.24670 =1 7202935 Py 2S 
ARH, , =a 9/3 OS -6.00543 1.2460 
ay =a a3 o9 =23 2.99202 2 oo 
ATI), =, 37943 —=51.,23575 Meares a Lear 
AHB,,. OF Oka 2 4 Os 0 1-4 4 Loe 4 
Constante —112260.12500 =112446-,18750 


iv bise Vel. SAN PDIEGO DRY SEASON FOG AND NO=FOG CLASSIFIcaA- 
TlONeEUNCTIONS FOR CATEGORIES 3 AND 4 


Ts 





COCTEVCLents. Of Jthe- 
VARIABLE FOG FUNCTION NO-FOG FUNCTION FINAL F VALUE 


































WET 1 ¢ (Removed at step 23) 

FP 16 O2274e7 =@.. 5.0 53341 Ons 54 

eB is -0.04832 =O 0 aaS ele Olay. 

ATIL, =e? Seek > OO =-25.54906 8.6549 

prBK oy OLS 13.5 | Ons 587) 326692 
A,TWB. (Removed at step 22) 

ASPD,. G67 513.356 LG. Sones 3 Deo 
A,DIR,, O- 01693 0.01492 325020 
A,Sl,, (Removed at step 19) 

TOtR G4 —-0.14458 -0.14257 1.9764 

ATWB,; -1.23822 -0.85878 4.6643 

ARH, ; ie B72 ies 202 Ped Ty 
A,TDB,, 20.99492 20.64488 6.13941 
ARH, -2.62901 -2.64525 2.293% 

12 16 wie to VO AZ See 10 5. 2060 

AsI,, -5.89893 =5 2525356 Bie 1 OG 
A,T1,, O72 90936 O37 25e0 BIS )S12h5) 

eat e : = eo e i 
ATI, Ses) (Si Alpayes 6.61983 PAP, 2), 8) 72) 
oT WB; Do oO 6 / oe Bape ea oi 0) 7 4239779 
rE, EL eso ki OM ON (Ens 2 SN3) 14.2448 
ok wee ae oiay 4.3569 
1 Lic to 7 70 90 
35116 -19.78639 -~19.68510 NEA ONS) 71 (5) 
ae e a ae ne eS Wee: 
\ 1.14815 Ie 2.05eu 
-O. =O. ~2 000 
AOHB 6 oa G70 13 67 0.01401 il O 

2°PDi HO. 54.596 10.48309 iS 206 

Sonmoabtant.—5/540.84766 -~-5/436.34766 


feo Ge VIL. SAN DIEGO WET SEASON FOG AND NO-FOG CLASSIFICA- 
PEON FUNCTIONS FOR CATEGORIES => AND =a 
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SEASON & CATEGORY HOGee ERCENTAGE NO-FrOG PERCE NIAC 





WEIGHTED AVERAGE 


Peooek VilIl. SAN DIEGO FOG AND NO-FOG DISCRIMINATION 
PERCENTAGES 
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CATEGORY NUMBERS.) CLASSIFICATION 


ROGGEKS 
NO={=2OGGERS 


BOGGERS 


NO-FOGGERS 





Pease IX. OAKLAND DRY SEASON FOG AND NO-FOG CASES 
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CATEGORY NUMBER  & CLASSIFICATION 


FOGGERS 
NO-F OGGERs 


FOGGERS 


NO-FOGGERS 





Peab GE Xx. OAKLAND WET SEASON FOG AND NO-FOG CASES 


a9 





Coefficients of the: 
VARIABLE FOG FUNCT TON NO-FOG FUNCTION FINAL F VALUE 


SPD 


A 
ail 


04 


OLE2G7 9 
=O" 010iz2 2 
1.05340 
-0.13446 
=O.22815 


(Removed 


O.07S02 


(Removed 


Oi 59e26 
O02 00227 
0.04460 
O20 0S 6 
Os63973 
O20e523 


0.06999 
=0 -OU0GG 
ie 6 0S 
O70 s 99 
On 0 G9su 
Step gis) 
OO. 22367 
step 16) 
O04 5:23 76 
—0 0.0 E05 
0.04094 
O.. 00256 
Or? S0eZ 
=0207332 


6.05 0n 
Sl ere: O10, 
4.5235 
3.24 
33 9Gec 


4,9993 


Pal ST) Fie, 
1.2546 
LO 2k 
1.6784 
Lev224 
SZ oe 


constant —2t5 309205 =1S 4462 





ABLE Xl. OAKLAND DRY SEASON FOG AND NO-FOG CLASSIFICATION 
PUNCTIONS FOR CATEGORIES I AND 2 
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CGoOoerinetents of the: 
VARIABLE FOGs"UNGCTION NO-FOG FUNCTION FINAL F VALUE 


"1G 


eG 


a O4 


ASI y, 


ee 16 


Bl 94-16 


A,WET,« 


16 
04 


A,DIR 
A, TDB 


ASI 


04-16 
ec 
A, THK, « 
ns 04 
eo 
ie 4 
oie! 16-04 


3 THK, 


a 


Constant 


tae ALT. 


3.34849 
= 3, 05567 
AZO 12s 0 
= Oe DOS 5 
=O 7 2a 

0.02788 
m0 IL Seis: 
«QO. P1382 


3B. 27436 
=O fo oS 
42.00005 
-9.84301 
= Oe 16 7 Sez 
~0.04618 
edsiey 0.2) 74 Ole: 
oO SILI OLS I5, 


(Removed at step 17) 


Sez Za 
=) aC bene, 
=0.09032 

Owls Oo 

O00 7e5 

On 02522 

LAO, 
= © 702065 

Jee J 


S27 517.31641 


OAKLAND WET SEASON FOG AND NO=-FOG CLASSIFICATION 
PUNGETONS FOR CA TEGORTES = eal Dar 


8.20044 
pO LONO! 2 2s 
10.7510) S) Tire) IL 

0.13412 

O20CS 97 

O.02869 

128298 
“QO. 0357 4 

9.94415 


= 2143). SOvene 
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TP Oa2 24 
A Sole eS: 
oie) IL I 
ae) 3), 
6. / 022 
(YS eZ 
220160 
4.4790 


225240 
2.6200 
Phra) IE 
1. 99a 
4./836 
3.00%7s 
Ze 5 oe 
HEPA 6) 10) IL 
Eel 70 








VARIABLE 


ADPD), 


AWBD,, 
aa 1G 
A,DIR 


Ne 
es 16204 
ee Old 
ae NG 
a 2 4-16 


ee 16 

5 WEG 

A,SLP Oy 
A a Ea 

ats 


5 vO4 


i 
ae 
ats 
1 ae 
ait 16 


ASPD), 


Constant 


Tees XLIL. 


GOocrttLLclenes of Ehe- 
FOG FUNCTION NO-FOG FUNCTION 


4.41206 Geo 9220 
mL) OLS. o)S 


O20 0352 


=o eo LOO 3 
OA LONG eyes, 
(Removed at step 12) 

=). 9130 One oe3 
=—1.24039 =o. O16 
1.40700 O27 73656 
(Removed Step 18) 
CL PS ciisy 23) 
he O09. oO 


0.30540 


-0.96936 

Coy 150 

-0.30712 
(Removed 


2.40190 


step 21) 
3.209862 
-0.65045 
2.62428 
-0,61391 
-0,.16509 
-0.60388 
H/o sos 954 


al Ougseey 2 
ie 5426 
seo 2 2 0 
0.44584 
-0.047/02 
=a Ges O08 


OAKLAND DRY SEASON FOG AND NO-FOG CLASSIFICATION 


BUNGCIIONS POR CATEGORIES. S anD. 4 


eZ 


FINAL F VALUE 


Bie BOG 
236275 
2230 6 


Cosa O Bone) 
Zo 9 
sso) 2) oe 


12.4420 


8.4864 
Zee 


SO 510 5 
ite 5 
e)ryts) LOE 
5.4850 
PLEATS ENS 
Sao 











Coertiacrents of the- 
VARIABLE FOG FUNCTION NO-FOG FUNCTION FINAL F VALUE 




































1RAi¢ 226776054 22.05 Gue 7.2474 
ities dl lS) Sys) (59 7 =i 43764 i 
35LP,, 2.84706 2.007 oe Seo oe 
ASPD, | = by ff 22) Jal ware SS) 01e:te: 10.4081 
3PIR,, 0.23442 OF.239 42 SS) SNS), 
A rr o4- 16 Ole S 3-4 1635 6.6105 
1H To -~0.01989 =O) OS) 1S 9.8690 
vr 16 64.52661 64.19128 ZOU GS 
37 BIG weeOS Oo 2 =2.48583 SB 01652 
oTDPi (Removed at step 23) 
3 THK), 0.02141 O02 208 3.4260 
A pW 04-16 Plog eA Sas us ass 7) 2.4149 
3F Soa =U es 2 7 1S i O61 
A,TT 04 -4.41379 -4.44339 POST 
A,SLP Oy a32.536093 ~32.94241 27 OG 
oFSi¢ =e Oo. 24606 —~ 104236405 eco? 
A FSi¢ oe eo 5 -57.47128 2. 9 96 
ADIR,, O258270 O23 5426) lito les 
AjPIRD A 16 0.21845 O22 cu 20623 
A,DPD,.« -0.53928 -0.87731 247606 
ARH "ué Wwe OGG 5 =I o2oo2 1.4863 
12 Bis ZO S 3:5 7 217 1) Gre 4.4619 
. = oso LOS0 a STS Aes) Mot GTS 
oSLP oO, 69.76450 GO Ge 762 Vee Oboe 





wonstant =—7Orl2 7/23/7891 -37993.56641 


OAKLAND WET SEASON FOG AND NO-FOG CLASSIFICATION 
HONECTIONS FOR CATEGORIES 3 AND ee 


TABLE X1V. 
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SAN & CATEGORY FOG PERCENTAGE NO-FOG PERCENTAGE 





WEIGHTED AVERAGE 


TABLE XV. OAKLAND FOG AND NO-FOG DISCRIMINATION PERCENTAGES 
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GATEGORY NUMBER & CLASSIFICATION 


FPOGGEKs 








NO-]FOGCGERS 


PmOCGE Rs 














NO-FOGGERS 


TABLE AVI. OULLEAYUTE DRY SEASON FOG AND NO=rOG CASES 
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CATEGORY NUMBER & CLASSIFICATION 


FOGGERS 
NO-FOGGERS 


FOGGERS 


NO[hOGGER. 





Ween XVIL. OQUILLLAYUTE WET SEASON FOG AND NO-FOG CASES 
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GOCEEMCI ONES (On. the. 
VARIABLE FOG FUNCTION NO-FOG FUNCTION FINAL F VALUE 


A,THK) (Removed at step 28) 
AWET,« (Removed at step 18) 
A, SPD, 6_o4 =12.59450 | -12.84844 
AjSI,, (Removed at step 23) 
git aa O. O02 7Ga 0.02604 34 4,0cal 
3° 416 23.46292 232108 7 220 oe 
3HB,.¢ OO ao7a9 O07 0565.0 345254 
ATI), =—954.500514 =e O orn 29.4205 
1P1Ro, -2.06415 -2.08957 21.8740 
A,RH, ¢ O02 7 54 =O 0 3 one> 4.8483 
A, SPD, 23,2594 72 23... 5 005 Det243 
A,DPD 4, oO TOw 263 55 10.0 Kew, 
W 7 
, BD a, (Removed at step 25) 
A,TBo, (Removed at step 29) 
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APPENDIX A 


A GENERALIZED FORTRAN PROGRAM FOR APPLICATION 
OF THESCLASSIFICATION BUNGEIONS 
VAR(I) = VARIABLE RAW VALUE 
N = NUMBER OF TERMS IN THE CLASSIFICATION EQUATION 
Crt) — FOG FUNCTION COEFFICIENT 
Chieti) — NO=-FOG FUNCTION COEFFICIENT 
DEP CG DISCRIMINANT SCORE 


Die NO-FOG DISCRIMINANT SCORE 


I! 


PPF POsrigntlOR PROBABILITY OF FOG 


PPNF = POSTERIOR PROBABILITY OF NO-FOG 
DIMENSION VAR(N),CN(N) ,CNF (N) 

pe =~ CONST 

DNG = CONST 

READ IN THE VARIABLE RAW VALUES 
READ(5,5) (VAR(I),I=1,N) 

PORMAT (8F10. 2) 

MUDBEDELY THE N VALUES OF THE VARIABLES BY THEIR 
MeSGecriATeD COBFFICIENTS AND SUM 

DO 10 I=1,N 

be =] DF + CF(1)*VAR(T) 

DNF = DNF + CNF(I)*VAR(I) 


SONTINUE 
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$516, 


TO AVOID TRYING TO EXPONENTIATE A NUMBER TOO LARGE 
BOeBE  KREPRESENTED BY THE IBM 360 GOMPUTER DO FHE 
FOLLOWING: 


DMIN AMIN1 (DF,DNF) 


DIFF = DF - DMIN 
tHe UEE.GH.-174.67) GO TO 26 

Deere = .0.0 

PF = EXP (DIFF) 

DIFF = DNF — DMIN 
PievMle bh. Gh.—1/4.6/) GO TO 30 

DIFF = 0.0 

PNF = EXP (DIFF) 

SP = PF + PNP 

NOW FORM THE POSTERIOR PROBABILITIES 
PPF = PF/SP 

PPNF = 1.0 - PPF 

NOW OUTPUT THE PROBABILITIES IN ANY DESIRED FORMAT 
STOP 


END 
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